Entanglement entropy of integer Quantum Hall states 



Ivan D. Rodriguez and German Sierra 
Instituto de Fisica Teorica, UAM-CSIC, Madrid, Spain 



o 
o 



I 

O 

o 



> 

00 
00 

(N 



oo 

o 



X 



We compute the entanglement entropy, in real space, of the ground state of the integer Quantum 
Hall states for three different domains embedded in the torus, the disk and the sphere. We establish 
the validity of the area law with a vanishing value of the topological entanglement entropy. The 
entropy per unit length of the perimeter depends on the filling fraction, but it is independent of the 
geometry. 



In recent years the notion of entanglement has become 
a new tool for analyzing the quantum states that arise in 
Condensed Matter systems This notion has brought 
a quantum information perspective to traditional prob- 
lems and techniques in the field, such as quantum phase 
transitions, numerical simulation methods, renormaliza- 
tion group, etc. A generic measure of entanglement is 
given by the von Neumann entropy Sa of the reduced 
density matrix of a part A of the total system. This 
quantity measures the amount of quantum entanglement 
of the subsystem A with its environment, usually de- 
noted as B. For finite systems, one has that Sa — Sb, 
so this quantity reflects a property shared by A and B. 
If the quantum state has a finite correlation length, then 
heuristic arguments implies that the entropy is propor- 
tional to the area of the common boundary between A 
and B. This statement is known as the area law, and it 
seems to be a universal property satisfied by the quan- 
tum states appearing in Nature (see Q for some rigor- 
ous results) . There are logarithmic violations of this law 
in critical one dimensional systems, and some higher di- 
mensional fermionic systems, but the former ones can be 
understood using Conformal Field Theory [3]. 

Particularly interesting are the two dimensional sys- 
tems with topological order, where the entropy law be- 
comes Sa = cL — 'y + 0{l/L), where L is the length of the 
boundary, c is a non universal constant and 7 is a quan- 
tity called topological entanglement entropy [1| . The ex- 
citations of these systems are anyons and it turns out 
that 7 is the logarithm of the total quantum dimension 
of these anyons. The paradigmatic system with (abelian) 
anyons is the Fractional Quantum Hall (FQH) state with 
filling fraction 1/m, for which 7=5 logm. The area, or 
rather, perimeter law of the FQH states has been the tar- 
get of several recent studies [1, @, 0] , in order to confirm 
its validity and to compute the value of 7 predicted in Q 
(see Q for the study of the entanglement entropy for 
particle partitioning). Reference [t^ uses Chern-Simons 
theory, finding the predicted value of 7, however the lin- 
ear behaviour of Sa, is not captured, due to the purely 
topological nature of this theory. There are numerical 
studies using the Laughlin wave function [5] and exact 
diagonalization for filling fractions i/ = l/3jl/5 and 
the = 5/2 Pfaffian state. The approaches of [1, @] use 
the orbital basis for the Landau levels. The close relation- 
ship of this basis to the spatial partitioning of the blocks 
leads to an area law of the form Sa = cy/TA — ^+0{l/lA), 



where Ia is the number of Landau orbitals in the block 
A. The numerical values of 7 computed in the spheri- 
cal geometry Q and the torus geometry Q agree, within 
some precision, with their theoretical values, despite of 
the fact that the systems analyzed are not very large. We 
remark that the previous form of the area law in the or- 
bital basis holds only in the case of fractional fillings. For 
integer fillings the orbital partitioning entropy is actually 
zero since the ground state is simply a product state in 
that basis. 

In this letter we address the problem of computing the 
entanglement entropy Sa directly in real space, for the 
Integer Quantum Hall states with > 1, in three dif- 
ferent domains: strips in the torus, annulus in the disk 
and casquettes in the sphere. The reason for choosing 
integer filling fractions is that the ground state is given 
by free ferniions, where standard techniques for comput- 
ing entanglement entropies are available Q . We find the 
area law Sa ~ Cj^L — 7, with 7 = in agreement with 
general arguments The non universal constant Ci, is 
computed analytically for v = 1 and numerically up to 
ly = 5. We also analyze the crossover from thin to large 
blocks, finding that the onset of the area law occurs when 
the width of the boundary is larger that a correlation 
length. The blocks, whose entropy we have computed, 
are adapted to the standard gauge choices used to ana- 
lyze the geometries of the torus, the disk and the sphere. 

Let us consider the Landau model for a particle in a 
torus of size x Ly. The one particle wave function in 
the lowest Landau level (LLL), in the gauge A — B{0, x), 
is (in units of the magnetic length i equal to one): 



1 



1/2 



pikyV „-{x-ky) /2 



(1) 



On the torus, the identification of the wave function along 
the y direction implies: 
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(2) 



The number of LLLs, no? is obtained imposing that the 
particle lives in the strip \x\ < Lxl2, which yields no = 
. This value also gives the total number of quantum 
fluxes through the box. The electron operator can be 
written as 



X^'V^ky higher LLs, 



(3) 
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where Ck^ is the fermionic destruction operator of the 
LLL labelled by ky. The extra term in ^ involves the 
remaining Landau levels, which are empty for filling frac- 
tion v — I. Later on, we shall take them into account 
when considering higher filling fractions u. The ground 
state for = 1 is given by: 
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(4) 



where |0) is the Fock vacuum. The two point fcrmion 
correlator in this state is. 

Cry = {^oH\x,y)i;{x',y')\^o)- (5) 
Using ^ and ^ one finds: 

Cr.r' =^0^, (x, 2/) (6) 

ky 

and plugging (P) 

Ky 

(7) 

The sum in ((7]) is over the no values of ky given in ([2]). 
In the limit L^, Ly ^ oo the correlator ^ becomes 

27r 

and it is short range with a correlation length propor- 
tional to the magnetic length £ = 1. 

We want to compute the entanglement entropy, St>, of 
the state $0: in the strip 

I?:-|<a;<|, 0<y<Ly. (9) 

This entropy is given by the formula 5p — Tri5^|$o)(<&o|, 
where Vc is the complement of V in the torus. The com- 
putation of St> is done in two steps @ . First one restricts 
the correlation matrix Cry, to the domain V, i.e. 

Cry ^ Cry, r,r'eV. (10) 

Next, one diagonalizes Cry, i.e. 

d^r' Cry f,n{r') = A,„/„,(r). (11) 

I 

The entropy Sx> is obtained by means of, 

Sv=Y.H{X,n), (12) 

m 

where i?(x) = — a: logx— (1— x) log(l— a;). The eigenvalue 
problem (jlip can be rather difficult for a generic domain 
V, however for the strip ([9]) this taks simplifies. The basic 
observation is that Cr.r' only depends on the difference 
y — y' , which suggests the ansatz 
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Plugging ((T3l) into (fill) , and taking the limit Ly ^ oo 
one gets 



where 



A — 



dx 



ix/2 



(14) 



(15) 



For a non vanishing eigenvalue Am , eq. (|14p fixes the func- 
tion gm{x), up to an overall factor. Plugging p4| into 
P^ . the constant A,^ drops, and one gets the eigenvalue 



dx 



-{x-tj,^f 



(16) 



On the other hand, if Am = , ea. p4)) yields Am = 0, 
which becomes a condition for the function g^. However, 
vanishing eigenvalues do not contribute to the entropy 
(fT^ . so the solution of Am = is not required. Recalling 
that the function ([13]) is defined on the domain ([9]), one 
obtains a quantization condition similar to eq.(l2|) 

tJ'm = ^—, , — ^ + l<m<y. (17) 



In fact, the eigenfunctions fm of Cr.r' coincide with the 
conjugate of the LLL eigenfunctions 0^ , under the iden- 
tification ky — /im- Moreover, eq. (|16p can be written as 
the norm of ([T]) over the domain ([9]), i.e. 



d^r |0^„(r)p. 



(18) 



which means that Am is the probability of finding the 
electron in the state ky = /im in the domain V. Integrat- 
ing yields 



Am = A(/im7^a;) — ^ 



Erf (/im + |) - Erf(/.m - |) 



(19) 

where Erf(x) is the error function. The function 
H{X{fi,lx)) is localized in the regions ^ Ix/"^, asso- 
ciated to the boundaries of T), where it can be approxi- 
mated as 



l_Erf(H-|) 



iErfc(|/i|-|), (20) 



where Erfc(a;) = 1 — Erf(x) is the complementary error 
function. In the limit Ly » 1, one can use eq. (fT7|) to 
write (fT^ as the integral, 



T f°° 

Sv = S{lx,Ly)^-^ dijiH{\{p.,lx)). 

J^oo 



(21) 



/m(r) = e-'^'^ygraix). 



(13) 



Furthermore, \i Ix » 1, the main contribution to ([2T|) 
comes from the values of ^ around ±lx/2, where one 



3 



S(l^,Ly) /Ly 

0.4 



. 3 

0.2 



. 1 



0.5 



1 . 5 



FIG. 1: Plot of S{lx,Ly) as a function of for Ly 
20, 30, 40. 



can use the approximation (j20p . Shifting the integration 
variable /i, one finally obtains 



where the constant Ctorus is given by 



(22) 



Ctorus = f ^ iJ(iErfcOi)) « 0.20329081. (23) 



Equation ([221) is the area law satisfied by Sv in the limit 
Ix >> 1. The absence of a constant term in (P^ . implies 
the vanishing of the topological entropy 7 [1]. 

Equations ((22)) and (j23p were derived under the asump- 
tion that >> 1, but using ea. (fT^ one can estimate 
the value of Ix above which the area law (|^^ starts to 
be valid. In fig. [U we plot S{lx, Ly) / Ly as a function 
of for various values of Ly. For Ix > Ix.c = 2 the 
entropy reaches a constant value 2ctorus given by equa- 
tion (P5)) , which agrees with the area law (|22p . The value 
of Ix^c is easy to understand, it means that the width 
of 2?, along the x direction, must be bigger than two 
magnetic lengths, which guarantees the existence of two 
boundaries with a width of one correlation length each. 
Fig. [ijalso shows that the entropy increases linearly with 
VT^ for Ix < lx,c, and reminds the result of reference [1], 
where the entropy of the ly — 1/3,1/5 FQH states was 
computed for the torus in the orbital basis. These au- 
thors found that the entropy of a strip of width £, varies 
as v^, where £ label the orbital angular momenta. Our 
computation and that of reference 6] differ both in the 
basis and the filling fraction, so the previous comparison 
has to be taken with care. 

Let us now consider the disk geometry in the Landau 
gauge where the eigenfunctions of the LLL are given by 



4>m{z) 



(27r2™m!)i/2 
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m = 0, 1, 



The two point fcrmion correlator is similar to eq. 



Cz.z' 
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(24) 



(25) 



We want to compute the entanglement entropy in the 
annulus of radii ri < r2 



V -.ri <\z\< r2. 



(26) 



The procedure follows closely the case of the torus. The 
eigenfunctions fm{z) of eq. (fTT|) . with non zero eigenvalues 
Xm, are given by fm{z) = cl)*^{z) (m = 0, 1, . . . ) and the 
eigenvalues are 



d?z |0m(2;)p, TO = 0, 1, 



(27) 
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Plugging ([M]) into (P7|) and performing the integral over 
the domain (l26l). one finds 



Am (n, 7-2) 



1 



r(m + i,^)-r(m + i,^) 



The entropy of the annulus is given by 

00 

S{ri,r2) = ^ H{\„,{ri,r2)), 



(28) 



(29) 



which, for large values of ri and r2, satisfies the area law 
5'(ri,r2) = 27r(ri -l-ra) Cdisk, (30) 

with Cdisk = Ctorus- Equation (j30p can be proved analyt- 
ically along the same lines as was done before. 

Another example which can be solved explicitely is 
that of an electron moving on a sphere of radius R un- 
der the infiuence of a radial magnetic field created by a 
monopole at the origin. In the gauge where the vector 
potential is given by A = HQ/eRcot (p, the wave func- 
tions of the LLLs are the monopole harmonics. 



Q,Q,m 



(2Q + 1 ( 2Q 
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(-1) 



Q-rn^Q+m 



(31) 



where u = cos(6'/2)e"*'^/^ v = sin(6'/2)e''^/^ with 9 and 
(j) the polar and azimuthal angles, and 2Q the total quan- 
tum fiux traversing the sphere. The two point fermion 
correlator is given by 



Cr,r' — ^ Yq, 



Ya 



>')■ (32) 



Making the change n = m + Q one can write ((32| as 

2Q 

i 
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= ^ E f ? ) (""T {vvT'^-. (33) 



ji=0 



We are interested in computing the entanglement entropy 
in the spherical segment (i.e. casquette) 



V:0a<e< 



(34) 
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bounded by the polar angles 6i and 02- The eigenfunc- 
tions of the correlator ([55)1 in the domain ([M]) , with non 
zero eigenvalue, are given hy fm — Yq q with 



A. 



de 



2-iT 



d(t>\yQ,Q,n- 



(35) 



Performing the integral one finds, 



^n\Oa^ Ob) 



X S(cos2(^),l + n,l-n + 2Q) 
B(cos2(^^),l + n,l-n + 2Q)j , (n = 0,l,. 



(36) 



,2g), 



The area law S*© « c^L — 7^ remains valid, with 7,^ = 0, 
and a value of c^, which depends on the filling fraction 
(see fig. [2|). 

In summary, we have derived in this letter the area 
law satisfied by the entanglement entropy of the integer 
Quantum Hall states with filling fraction, for different 
types of domains in the torus, the disk and the sphere. 
We have computed the non universal constant Cu of the 
area law as a function of v. The topological entangle- 
ment entropy vanishes, in agreement with the theoretical 



where B{z, n, m) is the incomplete beta function and 
^0(6) = a{b)n/ {2Q + 1). The entropy of the region ([34| is 
computed from 



2Q 



(37) 



The perimeter of the casquette of angle 9 is given by 
Pe — 2TTRsmd, where the radius is given by i? = y/Q as 
follows from computing the number of quantum fluxes. 
For large values of Q, the entropy ([57]) satisfies the area 
law 




FIG. 2: The points denote the values of the constant Ci, in 
the area law for integer filling fractions v — 1,...,5. The 
continuous line is a forth order polynomial fit. 



S{0a,0b) = {P{0a) + P{Ob)) Csphcrc; (38) 

with Cgphere = Ctorus — Cdisk, SO that the three geometries 
yield the same entropy per unit length of the perimeter. 

The previous results can be easily generalized for inte- 
ger filling fractions v > 1. The correlation matrix Cr,r' 
is given by 

Cry^ E E'^".™W'^«,™(r'), (39) 

n— 0,1/— 1 m 

where (t>n,m(j) is the wave function of the state m in the 
n*^ Landau level. The eigenfunctions of C^.r' are linear 
combinations of 0* ^(r) with n — 0,1, ... ,1^ — 1, and 
their eigenvalues are those of the v x v matrix, 

Arn{n,n')^ / d^r (/)* „(r)0„/.„(r), n, n' = 0, . . . , 

(40) 

The entanglement entropy is computed using ea. (|12p . 
where the summation runs over all the eigenvalues of A^. 



results [J, |7|. For = 1, we have found a simple inter- 
pretation of the area law. In this case the entanglement 
entropy is given by the sum over the LLL states, of the 
Shanon entropies associated to finding an electron or a 
hole in the domain. The area law arises from the con- 
tribution of the LLLs inside a correlation length of the 
boundary of the domain. Our method allows the compu- 
tation of the entanglement entropy for more complicated 
domains. Of special interest are those with curvature 
singularities wether one may expect deviations from the 
area law. 
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